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Beschreibung 

Verfahren und Anordnung zur Ermittlung von Par ame tern eines 
technischen Systems 

Die Erfindung betrifft ein Verfahren sowie einer Anordnung zur 
Ermittlung von Parametern eines technischen Systems 

Bei einer Vielkanal-Obertragung und einem Vielkanal-Empf ang von 
Signalen tritt h&ufig eine Interferenz beispielsweise zwischen 
den Signalen/Bildern auf . Ein typisches Beispiel ist dabei eine 
Mischung eines Sprachsignals mit Rauschen, was bei der 
Telekoramunikation und bei Videokonferenzen ein grofies Problem 
darstellen kann. Die vorliegende Erfindung betrifft daher das 
Gebiet der Signaltrennung, um beispielsweise ein urspriingliches 
Sprachsignal wiederzugewinnen, 

Typische bekannte Techniken zur Trennung der Quellensignale 
ausgehend von Mischsignalen basieren auf einer zeitlichen 
Mittelung oder einer Filterung der Signale. Dies hat indessen 
Nachteile beztiglich des Rechenauf wands zur Folge. 

Es sind auch Verfahren auf der Grundlage der sogenannten Blind 
Channel Equalization (Signalentzerrung ohne Vorkenntnisse liber 
den Obertragungskanal) bekannt, aber diese Verfahren benStigen 
immer eine gewisse Kenntnis iiber die Quellsignale, wie 
beispielsweise eine Kenntnis iiber ihre statistische Verteilung. 

Das Problem der Signaltrennung tritt beispielsweise auch auf, 
wenn zwei Sprecher in zwei entfernt von ihnen stehende 
Mikrofone sprechen, so dafi jedes Mikrofon ein Gemisch der von 
den zwei Sprechern gesprochenen Signale aufnimmt. Somit besteht 
das Problem, das Signalgemisch, also eine Menge Uberlagerter 
Eingangssignale wieder zu trennen* Aus L.Molgedey, H.G. 
Schuster, ^Separation of a Mixture of Independent Signals using 
Time-Delayed Correlations", Phys. Ref. Lett. 72, 3634 (1994) 
ist dabei das folgende Verfahren bekannt: Das Problem, n 
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uberlagerte und korrelierte Quellensignale (Eingangssignale) zu 
trennen und gleichzeitige die Mischungskoeffizienten der 
Quellenstarken zu bestimmen, laflt sich auf ein Eigenwertproblem 
reduzieren, bei dem simultan zwei symmetrische n x n Matrizen 
diagonalisiert werden mtlssen. Die Matrixelemente sind meBbare 
zeitverzSgerte Korrelationsfunktionen. Die Losung des 
Eigenwertproblems kann durch ein neuronales Netz erfolgen, 
wobei die Lernregel fttr die lateralen inhibitorischen 
Wechselwirkungen zwischen den Neuronen durch eine Liapunov- 
Funktion bestimmt werden kann, deren Minima die (entarteten) 
Losungen des Problems lief em. 

Dieses Verfahren wurde auch bereits auf akustische 
Eingangssignale angewandt (siehe F. Ehlers, H. G. Schuster, 
„Blind Separation of convolutive mixtures and an application in 
automatic speech recognition" IEEE Trans. Signal Proc. (1997) . 

Aus der DE 195 31 388 CI ist ein Signaltrennungsverfahren und 
eine Signaltrennungseinrichtung fur nichtlineare Mischungen 
unbekannter Signale (Blind Channel) bekannt, das schematisch in 
Fig. 3 dargestellt ist. 

Dieses deutsche Patent behandelt die Separierung eines 
Signalgemisches, bestehend aus der nichtlinearen Oberlagerung 
von M unbekannten Quellsignalen XI, X2, wobei N (N £ M) 
unterschiedliche Mischungen der M Quellsignale XI, X2 inklusive 
eines eventuellen Storsignals einer Signalauswerteeinrichtung 
zugeftihrt werden, die durch eine statistische Analyse der 
Signale die nichtlinearen Obertragungsfaktoren bestimmt und mit 
diesen errechneten Faktoren die Mischung rUckgangig macht, so 
dafi die N Ausgange der Signaltrennungseinrichtung moglichst 
nSherungsweise die M Quellsignale ohne Oberlagerungen 
enthalten. Dadurch wird eine Behandlung nichtlinearer Gemische 
mOglich, wobei nichtlinear bedeutet, dafi die Quellsignale XI, 
X2 durch ein unbekanntes nichtlineares System G gemischt 
werden. Das unbekannte System G wird durch eine sogenannte 
Volterra-Reihe beschrieben, und die Signaltrennungseinrichtung 
G-l bestimmt die Koef fizienten der Volterra-Reihe. Mit deren 
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Kenntnis 1st eine Entmischung des Signalgemischs mogl^ch. 
Auflerdem ktinnen die Koeffizienten zu weiteren Analysen zur 
Orts- Oder Geschwindigkeitsbestimmung der Signalquellen benutzt 
werden . 

Das aus dieser Druckschrift bekannte Verfahren besteht dabei im 
wesentlichen aus zwei Schritten: 

• Erstens werden die durch den wahlbaren Grad der 
Nichtlinearitat bei der Mischung eindeutig bestimmten 
nichtlinearen Gleichungen ftir ein gleitendes Zeitfenster 
gelost und die Losungen werden tiber die Zeit gemittelt. Diese 
zeitliche Mittelung stellt einen Hauptnachteil dieser 
bekannten Technik dar, da sie den Rechenaufwand und 
gleichzeitig die Zeitdauer fiir die Berechnung erhoht. 

• Zweitens wird aus einer gentigend grofien Anzahl von 
unterschiedlichen Kumulanten der geschatzten Ausgangssignale 
gebildetes Potential minimiert, wobei die zur Berechnung des 
Potentials nStigen Werte aus einem gleitenden Zeitfenster 
einer wahlbaren LSnge stainmen. Dabei ist vorausgesetzt, dafl 
sich das Mischungssystem so langsam Sndert, dafl diese 
Anderung bei der Berechnung der gesuchten Mischungsfaktoren 
vernachlassigt werden kann. Gemafi diesem deutschen Patent 
wird zur Ausftthrung des zweiten genannten Schrittes eine 
Kostenfunktion konstruiert, die minimiert wird. Wenn das 
globale Minimum erreicht ist, hat man die optimalen Werte, in 
diesem Fall der Obertragungsfaktoren, gefunden. 

Hinsichtlich des zeit lichen Auf wands sowie des Rechenauf wands 
ist das in der DE 195 31 388 CI beschriebene Verfahren aufgrund 
der Ausfiihrung einer zeitlichen Mittelung am AbschluB des 
ersten, oben genannten Verf ahrensschrittes nachteilig. 

Die vorliegende Erfindung hat daher zur Aufgabe, ein Verfahren 
und eine Anordnung bereitzustellen, die die Trennung 
ttberlagerter, statistisch voneinander unabh&ngiger akustischer 
Signale mit verringertem Rechenaufwand ermOglicht. 
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Diese Aufgabe wird durch das Verfahren sowie durch die 
Anordnung mit den Merkmalen gemSJJ den unabhangigen Ansprttchen 
gel5st. 

Bei einem Verfahren zur Ermittlung von Parametern eines 
technischen Systems, mit dem Ausgangssignale aus einer Menge 
tiberlagerter, statistisch voneinander unabhSngiger 
Eingangssignale ermittelt werden konnen, werden die Parameter 
derart ermittelt, dafc die statistische Unabhangigkeit der 
Ausgangssignale maximiert wird. 

Eine Anordnung zur Ermittlung von Parametern eines technischen 
Systems, mit dem Ausgangssignale aus einer Menge iiber lager ter, 
statistisch voneinander unabhangiger Eingangssignale ermittelt 
werden konnen, weist einen Prozessor auf, der derart 
eingerichtet ist, dali die Parameter derart ermittelt werden 
kSnnen, da6 die statistische Unabhangigkeit der Ausgangssignale 
maximiert wird. 

Vorteilhafte Weiterbildungen der Erfindung ergeben sich aus den 
abhangigen Anspruchen. 

Die Parameter werden bevorzugt in einem iterativen Verfahren 
ermittelt. 

In einer weiteren Ausgestaltung sind die Parameter Elemente 
einer Entmischmatrix, mit der die Menge der uberlagerten 
Eingangssignale multipliziert oder auch gefaltet wird, wodurch 
die Ausgangssignale gebildet werden* 

Die Optimierung der Parameter der Entmischmatrix wird bevorzugt 
durch die folgenden Schritte erhalten: 

- Wiederholung einer zeitverzSgerten Dekorrelationsberechnung 
zur Ermittlung der Eigenwerte der Entmischmatrix, 
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- Ermittlung der Eigenwerte der Entmischmatrix, ftir die 
Kreuzkorrelationen einen minimalen Wert annehmen, und 

- Ausftihrung einer Kumulantenndnimierung, wobei als Startwerte 
ftir die Kumulantenminimierung die im vorherigen Schritt 
ermittelten Eigenwerte verwendet werden. 

Die Kumulantenminimierung kann beispielsweise durch Trainieren 
eines neuronalen Netzes, aber auch jegliche andere bekannte 
Minimierungstechnik, wie beispielsweise Gradientenabstieg oder 
Monte-Carlo-Simulationen verwendet werden. 

In einer Weiterbildung wird bei der Optimierung der Parameter 
der Entmischmatrix wenigstens ein Diagonalpararaeter der 
Entmischmatrix auf einen vorgegebenen Wert gesetzt, womit , um 
eine Stabilitat des Minimierungsvorgangs hin zu einem globalen 
Minimum zu gewahrleisten. 

Die Entmischmatrix wird bevorzugt auf eine finite Impulsantwort 
begrenzt, d,h. es wird ein FIR-Filter (Finite impulse Response) 
eingesetzt zur Bildung der einzelnen Komponenten der 
Entmischmatrix. Der FIR-Filter kann sowohl ein kausales FIR- 
Filter Oder auch ein nicht-kausales FIR-Filter sein. 

Ferner wird die Entmischmatrix wShrend der 
Kumulantenminimierung bevorzugt durch Projektion in einen 
Einheitskreis stabilisiert . 

Die Weiterblldungen gelt en sowohl ftir das Verfahren als auch 
ftir die Anordnung, wobei bei der Anordnung jeweils der 
Prozessor derart eingerichtet ist, daB der entsprechende 
Verfahrensschritt durchftihrbar ist oder durchgeftihrt wird. 

Die Erfindung sowie deren Weiterbildungen kdnnen vorteilhaft 
eingesetzt werden zur Trennung tiberlagerter, statistisch 
voneinander unabh^ngiger Eingangssignale, insbesondere 
akustischer Eingangssignale. 
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Das Verfahren sowie die Anordnung sind far eine beliebige 
Anzahl von Eingangssignalen anwendbar. 

Weitere Vorteile, Merkmale und Eigenschaften der vorliegenden 
Erfindung werden nunmehr bezugnehmend auf die beiliegenden 
Figuren der Zeichnungen anhand eines Ausfahrungsbeispiels naher 
erl&utert. 

Fig. 1 zeigt die Anwendung eines Systems zur Trennung 
Oberlagerter, statistisch voneinander unabh&ngiger 
akustischer Signale gemafi dem Ausfuhrungsbeispiel, 

Fig. 2 zeigt eine symbolische Darstellung des Systems von 
Fig. 1, und 

Fig. 3 zeigt eine aus dem Stand der Technik <DE 195 31 388 
CI) bekannte Signaltrennungseinrichtung far nichtlineare 
Mischungen unbekannter Signale. 

Zur Wiedergewinnung des urspriinglichen Sprachsignals wird aus 
einer Mischung von Signalen die statistische Unabhangigkeit 
zwischen den Quellensignalen (ursprangliches Sprachsignal und 
Rauschen), im weiteren auch als Eingangssignale bezeichnet, 
ausgenutzt und der inverse Vorgang des dynamischen Systems, der 
die Mischung der Signale ergeben hat, wird im wesentlichen 
naherungsweise trainiert (gelernt) . Zwei verschiedene 
Mischungen des Sprachsignals bzw. des Rauschsignals werden 
beispielsweise durch zwei Mikrofone 1, 2 {vgl. Fig.l) erhalten, 
die voneinander beabstandet sind und/oder in entgegengesetzten 
Richtungen ausgerichtet sind. Bei dem Verfahren wird die 
sogenannte zeitverzdgerte Dekorrelationstechnik (TDD, time 
delayed decorrelation) verwendet, urn die Lernphase zu 
initiieren, d.h. urn Startwerte fUr die Lernphase zu ermitteln 
und vorzugeben, wodurch der Berechnungsaufwand far eine im 
weiteren beschriebenen Kumulantenminimierung verringert werden 
kann und die Gefahr lokaler Minima verringert werden kann. 
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Fig.l zeigt zwei Mikrofone 1, 2, die ein erstes Eingangssignal 
Zl(t) und ein zweites Eingangssignal Z2(t) aufnehmen. Diese 
Eingangssignale Zl(t) und Z2(t) kfinnen untereinander wiederuin 
jeweils mit Rauschen vermischt sein, was durch eine Mischraatrix 
S (siehe Bezugszeichen 3) symbolisch in Fig.l dargestellt ist. 
Nach dem Empfang bzw. der Obertragung wird eine Menge XI (t) und 
X2(t) iiberlagerter, statistisch voneinander unabhangiger 
Eingangssignale Zl (t) und Z2(t) erhalten. Diese Signale werden 
in eine Berechnungseinheit 4 eingegeben, in der im wesentlichen 
zwei Schritte ausgefiihrt werden, die symbolisch durch eine 
Berechnungseinheit B (Bezugszeichen 6) far den ersten Schritt 
sowie ein neuronales Netzwerk 5 filr den zweiten Schritt 
dargestellt ist. 

Durch die Berechnungseinheit 4 werden zwei Ausgangssignale 
Yl(t) bzw. Y2(t) ermittelt, die bei optimaler Einstelllung der 
Parameter in der Berechnungseinheit 4 naherungsweise gleich den 
Eingangssignalen Zl(t) bzw. Z2(t) sind. Mit anderen Worten, bei 
optimaler Einstellung der Parameter der verwendeten Matrizen in 
der Berechnungseinheit 4 erfolgt durch diese Berechnungseinheit 
4 im wesentlichen der inverse Vorgang zu dem dynamischen 
Mischvorgang, der durch die Matrix S (Bezugszeichen 3) 
symbolisch dargestellt ist. Das Ausfuhrungsbeispiel beschaftigt 
sich mit dem Optimierungsvorgang der Einstellung der Parameter 
der Entmischmatrix. 

Es werden die Parameter der Matrizen in der Berechnungseinheit 
4 derart optimiert, dafi die statist ische UnabhSngigkeit 
zwischen den durch den Matrizenvorgang in der 
Berechnungseinheit 4 gewonnenen Ausgangssignale Yl (t) , Y2(t) 
maximiert wird. Zu diesem Zweck werden die Ausgangssignale 
Yl(t) bzw. Y2(t) zu der Berechnungseinheit 4 zurUckgefUhrt (s. 
RttckfUhrschleifen 7 bzw. 8). Durch ein iteratives Verfahren 
wird ermittelt, ob sich die statistische UnabhSngigkeit der 
Ausgangssignale Yl(t) bzw, Y2(t) im Vergleich zu dem vorherigen 
Schritt der Iteration erhoht hat (und somit die Iteration die 
„richtige* Richtung in Richtung des globalen Minimums einer im 
weiteren beschriebenen Kostenfunktion einnimmt) oder nicht. 
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Pig. 2 zeigt eine mathematische Darstellung des Schemas von 
Fig.l, wobei der Mischvorgang 3 durch eine Matrix S{q) 
mathematisch beschrieben werden kann und der Entmischvorgang, 
der durch die Berechnungseinheit 4 erfolgen soil, durch eine 
Entmischmatrix M(q) symbolisiert wird. 

In Fig. 2 ist somit das Problem der Trennung einer sogenannten 
Multi -Channel Blind Source (Vielfach-Kanal-Quelle ohne a- 
prior i-Kenntnis) in zwei Dimensionen dargestellt. Dabei ist 
angenommen, daB das Mischsystem S (q) , wobei q f tir eine 
EinheitsverzCgerung steht, stabil ist und gleichzeitig auch 
eine stabile Invertierung aufweist, d.h. daB es ein 
Minimalphasensystem ist. Daruber hinaus ist angenommen, daB die 
Eingangssignale Zl(t) und Z2(t) (beispielsweise ein Sprach- 
bzw, ein Rauschsignal) statistisch voneinander unabhSngig sind 
und nicht gaussfdrmig verteilt sind. Die Menge XI (t) und X2(t) 
der Uberlagerten Eingangssignale Zl{t) und Z2(t) sind 
Eingangssignale in ein Entmischsystem mit einer Entmischmatrix 
M(q), deren Parameter (Matrixelemente) auf eine Maximierung der 
statistischen Unabhangigkeit zwischen den Ausgangssignalen 
Yl(t) und Y2{t) trainiert werden, Unter „Trainieren w ist dabei 
der gut bekannte Lernvorgang beispielsweise eines neuronalen 
Netzes bekannt, das als ein Beispiel fttr eine Technik genannt 
sein soil, die statistische UnabhSngigkeit zu maximieren. Dies 
erfolgt durch Minimierung einer im weiteren beschriebenen 
Kostenfunktion J(M). 

Es wird eine Kumulanten-Kostenfunktion gebildet, die die 
Diagonalkumulantenelementen der Kumulantenordnung 2-4 
minimiert: 

4 2 
Dcum « J(M) = £ £ |c (i) nondiagj 

1=1 nondiag 

Folgende Gesichtspunkte der dynamischen Mischung durch die 
Mischmatrix S(q) sind dabei zu berxicksichtigen: 
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Stabilitat des Entmischsystems: 
Dies wird erreicht, wenn M(q) auf eine finite 
Impulsantantwort (FIR-Filter) beschrSnkt wird. Die 
Stabilitat des FIR-Systems M(q) kann darUber hinaus auch 
dadurch erhalten werden, dafi wShrend der Lernphase eine 
Projektion in den Einheitskreis erfolgt Eine 
m6glicherweise vorliegende Nichtkausalitat der 
Invertierung von S(q) kann durch eine geeignete 
Zeitverschiebung (Verz6gerung) des Eingangssignals X(t) 
kompensiert werden. 

Eindeutigkeit der getrennten Signale Y(t): 
Ftlr den Fall statischer Mischungen werden die 
ursprtinglichen Quellsignale durch eine Skalierung 
wiedergewonnen. Far den Fall einer dynamischen Entmischung 
ist die Gefahr einer Nichteindeutigkeit der getrennten 
Signale Y(t) sogar noch grdfier. Es ist of fensichtlich, dafi 
fur den Fall, dafi Yl(t) und Y2(t) statistisch voneinander 
unabhangig sind, auch jegliche linear gefilterte 
Modifikation dieser Signale immer noch statistisch 
unabhMngig sind. Daher ist eine zusatzliche Information 
notwendig, um die inharente Nichteindeutigkeit des 
Problems zu verringern. 

Gaussverformung der Daten: 

Algorithmen auf Kumulantenbasis flir eine statische Blind 
Source-Trennung minimieren bzw. eliminieren in effektiver 
Weise Diagonalkumulanten haherer Ordnung entsprechend den 
Ausgangssignalen Y(t). Andererseits* ftihrt eine lineare 
Filterung zu einer Gaussverteilungs-Verformung der Daten, 
bei denen die Kumulanten hoherer Ordnung in Richtung 0 
gehen. Dies kann daher zu Randltfsungen fUhren, in denen 
die Kostenfunktion ein lokales Minimum erreicht, wobei die 
erwiinschte eigentliche Trennung (globales Minimum) nicht 
erfolgt. Um diesen ungewtinschten Fall zu vermeiden, wird 
die Struktur der Entmisch-Transferfunktion 
(Entmischmatrix) M(q) einigen Beschrankungen unterworfen. 
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Urn die obengenannten Probleme zu umgehen, wird ein Ansatz 
gewShlt, daB wenigstens eines (Oder auch alle) der 
Diagonalelemente auf den Einheitswert gesetzt werden: 

Mll(q) = 1 
und 

M22{q) « 1. 

Diese Annahme ist exakt, wenn die Mischelemente Sll (q) und/oder 
S22(q} ebenfalls einen Einheitswert „V y aufweisen. Sonst sei 
angenommen, daB Sll(q) und/oder S22 (q) stabile Invertierungen 
aufweisen, was die Skalierung der Diagonalelemente von M(q) auf 
den Einheitswert erlaubt. Dieser Ansatz verringert wesentlich 
die Nichteindeutigkeit einer LSsung und vermeidet daruber 
hinaus in effektiver Weise die Gefahr einer ubermafiigen 
Gaussverteilungs-Verformung der Ausgangssignale. Auch wenn auf 
den ersten Blick die BeschrSnkung der Diagonalelemente von 
M(q), wie oben ausgefiihrt, sehr restriktiv wirkt, ist diese 
Annahme in der praktischen Anwendung in der Regel erftillt. Ein 
typisches Beispiel ist die Entrauschung von Sprachsignalen auf 
der Grundlage einer Aufzeichnung mit zwei Mikrofonen, wobei die 
Mikrofone voneinander raumlich getrennt sind oder ein Mikrofon 
in Richtung des Sprechers gerichtet ist, wShrend das andere 
Mikrofon in der umgekehrten Richtung gerichtet ist, so daB das 
zweite, von dem Sprecher abgewandte Signal im wesentlichen nur 
ein Rauschsignal aufnimmt. 

Der Kumulantenahsatz basiert auf einer direkten Ermittlung der 
Diagonalkumulanten, wie sie in dem eingangs genannten Artikel 
von F. Ehlers und H. Schuster, „Blind Separation of convolutive 
mixtures and an application in automatic speech recognition" 
IEEE Trans. Signal Proc. (1997), ausgefahrt ist, was aber 
inherent den Nachteil aufweist, daB die numerische LSsung sehr 
aufwendig ist. Daher wird eine geeignete Initialisierung dieses 
Minimi erungsverfahr ens verwendet. Urn Startwerte fur das 
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Minimierungsverfahren zu ermitteln, wird geraSB der vorliegenden 
Erf indung die Technik der zeitverzogerten Dekorrelation (TDD) 
zur gleichzeitigen Dekorrelierung zweier verschiedener 
ZeitverzSgerungen angewendet, wobei diese TDD-Technik sich auf 
ein geeignetes Matrix-Eigenwertproblem zuriickfahren laBt. Wie 
bereits gesagt, wird diese TDD-Technik gemSB der vorliegenden 
Erf indung zur Initiierung des Diagonal (Kreuzkorrelations) - 
Kumulantenminimierungsproblems verwendet • 

Zusammenfassend laflt sich das Verfahren in die folgenden beiden 
Schritte unterteilen: 

1. Wiederholung des TDD-Verfahrens auf Grundlage des 
Eigenwertproblems im Frequenzbereich far verschiedene 
Verz6gerungspaare und Ermittlung der Lbsung, bei der die 
Kreuzkorrelationsterme einen minimalen Wert aufweisen. 

2. Initiierung (Start) der Diagonal-Kumulantenminimierung auf 
Grundlage der im oben gen'annten Schritt ermittelten 
Startwerte { FIR- Parameter ) . 

Im folgenden sollen noch einige Haupteigenschaften und Vorteile 
noch einmal zusammengefaBt werden: 

- Es ist kein a-priori-Wissen der Signaleigenschaften 
notwendig, mit der Ausnahme, daB eine statistische 
UnabhSngigkeit gefordert wird. 

- Die Stabilitat des dynamischen Entmischsys terns wird durch 
die Modulierung seiner Komponenten als FIR-Filter 
gewShr leistet . 

- Eine ttbermaMge Gaussverteilung-Verformung wird durch den 
Ansatz vermieden, daB wenigstens eines der Elemente der 
Mischtransfer-Funktionsmatrix (Entmischmatrix) auf den 
Einheitswert gesetzt wird bzw, auf den Einheitswert 
skaliert werden kann, und 
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- da der Kumulanten-Minimierungsschritt (Schritt 2) einen 
hohen Rechenaufwand erfordert, wird der Lernalgorithmus 
beispielsweise eines neuronalen Netzes durch das TDD- 
Verfahren initialisiert . 
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Im weiteren ist ein Prograinm in Matlab, Version 4 oder Version 
angegeben, mit dem das oben beschriebene AusfUhrungsbeispiel 
auf einem Rechner realisiert werden kann: 

function [ cos t , out 1 , out 2 ] =cumulant_costFIRa2 (par, input , pi , p 1 1 , 
p2, p22,a3,a4) ; 

% [ cos t , out 1 , out 2 ] as cumulant_costFIRa2 (par , input , pi , pi 1 , p2 , p22 , 
a3 f a4) ; 

% cumulant cost 

% FIR representation used 

% filter function used in both directions (non_causal) 

[np,mp]=size (par) ; 

firl=par (l:pl) ; 

f i r 1 l«par ( 1+pl : pl+pl 1 ) ; 

fir2=par (1+pl+pll :pl+pll+p2) ; 

fir22=par (pl+pll+p2+l:mp) ; 

den-1; %FIR only 

out 1= [ input ( : ; 1 ) -filter ( f ir 1 , den, input ( : , 2 ) ) -f lipud ( filter 
([0 firll ], [den],flipud(input(:,2)}))]; 

%/std (input (1, : ) ) ; %dlsim 

%filter 

out2= [input ( : , 2) -filter (fir2, den, input (:,!)) -f lipud (filter 
([0 fir22 ] f [den] r flipud{input{: # l))))];%/std(input(:,2)J; 

%dlsim %filter 

out=[outl out2]; 

%outl ss outl/std(outl) ; % this scaling was not needed in examples 
in SIP98 paper 
%out2=out2/std(out2) ; 

Ld=0; % number of delays in calculating the cross-correlation 

cost3«0; 

cost4=0; 
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COStALLl=[3 ; 
costALL2=[J; 
ol2=outl.*out2; 
cost2=roean(ol2) A 2; 

oll2-outl.*outl.*out2; 
ol22=out 1 . *out2 . *out2 ; 

if a3 =1 

cost3=[mean (ol22) ] A 2+ [mean (ol22) ] A 2; 

end 

if a4 ==1 

cost4= [mean (oll2 . *out 1 ) -3*mean (out 1 . A 2 ) *mean (ol2 ) ] 
[mean( (outl. A 2) .*(out2. A 2) ) -2*mean (ol2) A 2- 
mean(outl. A 2)*mean(out2. A 2)] A 2+. . . 
[mean (ol22 . *out2 ) -3*mean (out2 . A 2 ) *mean ( ol2 ) ] A 2 ; 

end 

%cum4a=[cum4x(outl,outl,outl,outl) ] A 2; 
%cum4b=[cum4x(out2, out 2, out 2, out 2) ] A 2; 



cost=cost2+a3*cost3+a4*cost4; 



%-cum4a-cum4b; 
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Patentanspruche 

1. Verfahren zur Ermittlung von Parametern eines technischen 
Systems, mit dem Ausgangssignale aus einer Menge ttberlagerter, 
statistisch voneinander unabhSngiger Eingangssignale ermittelt 
werden kSnnen, bei dem die Parameter derart ermittelt werden, 
dafi die statistische UnabhSngigkeit der Ausgangssignale 
maximiert wird. 

2. Verfahren nach Anspruch 1, 

bei dem die Parameter in einem iterativen Verfahren ermittelt 
werden. 

3- Verfahren nach Anspruch 1 oder 2, 

bei dem die Parameter Elemente einer Entmischmatrix sind, mit 
der die Menge der ttberlagerten Eingangssignale multipliziert 
wird, wodurch die Ausgangssignale gebildet werden. 

4. Verfahren nach Anspruch 3, 

bei dem die Optimierung der Parameter der Entmischmatrix durch 
die folgenden Schritte erhalten wird: 

- Wiederhclung einer zeitverzdgerten Dekorrelationsberechnung 
(6) zur Ermittlung der Eigenwerte der Entmischmatrix, 

- Ermittlung der Eigenwerte der Entmischmatrix, fur die 
Kreuzkorrelationen einen minimalen Wert annehmen, und 

- Ausftthrung einer Kumulantenminimierung (5), wobei als 
Startwerte ftir die Kumulantenminimierung die im vorherigen 
Schritt ermittelten Eigenwerte verwender werden. 

5. Verfahren nach Anspruch 4, 

bei dem die Kumulantenminimierung durch Trainieren eines 
neuronalen Netzes (5) erfolgt. 

6. Verfahren nach einem der AnsprUche 3 bis 5, 



GR 98 P 2958 



16 

bei dem bei der Optimierung der Parameter der Entmischmatrix 
wenigstens ein Diagonalparameter der Entmischmatrix auf einen 
vorgegebenen Wert gesetzt wird. 

7. Verfahren nach einem der AnsprUche 3 bis 6, 

bei dem die Entmischmatrix auf eine finite Impulsantwort 
begrenzt wird. 

8. Verfahren nach einem der Anspriiche 3 bis 7, 

bei dem die Entmischmatrix wahrend der Kumulantenminimierung 
(5) durch Projektion in einen Einheitskreis stabilisiert wird. 

9. verfahren nach einem der Anspriiche 1 bis 8, eingesetzt zur 
Trennung uberlagerter, statist isch voneinander unabhangiger 
Eingangssignale . 

10. Verfahren nach einem der AnsprUche 1 bis 8, eingesetzt zur 
Trennung uberlagerter, statistisch voneinander unabhangiger 
akustischer Eingangssignale. 

11. Anordnung zur Ermittlung von Parametern eines technischen 
Systems, mit dem Ausgangssignale aus einer Menge Uberlagerter, 
statistisch voneinander unabhangiger Eingangssignale ermittelt 
werden kbnnen, mit einem Prozessor, der derart eingerichtet 
ist, daB die Parameter derart ermittelt werden konnen, daB die 
statistische Unabhangigkeit der Ausgangssignale maximiert wird. 

12. Anordnung nach Anspruch 11, 

bei der der Prozessor derart eingerichtet ist, daB die 
Parameter in einem iterativen Verfahren ermittelt werden. 

13. Anordnung nach Anspruch 11 oder 12, 

bei der der Prozessor derart eingerichtet ist, daB die 
Parameter Elemente einer Entmischmatrix sind, mit der die Menge 
der uberlagerten Eingangssignale multipliziert wird, wodurch 
die Ausgangssignale gebildet werden. 
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14. Anordnung nach Anspruch 13/ 

bei der der Prozessor derart eingerichtet ist, daB die 
Optimierung der Parameter der Entmischmatrix durch die 
folgenden Schritte erhalten wird: 

- Wiederholung einer zeitverzogerten Dekorrelationsberechnung 
(6) zur Ermittlung der Eigenwerte der Entmischmatrix, 

- Ermittlung der Eigenwerte der Entmischmatrix, fttr die 
Kreuzkorrelationen einen minimalen Wert annehmen, und 

- Ausfuhrung einer Kumulantenminimierung (5), wobei als 
Startwerte far die Kumulantenminimierung die im vorherigen 
Schritt ermittelten Eigenwerte verwendet werden. 

15. Anordnung nach Anspruch 14, 

bei der der Prozessor derart eingerichtet ist, daB die 
Kumulantenminimierung durch Trainieren eines neuronalen Netzes 
(5) erfolgt. 

16. Anordnung nach einem der AnsprUche 13 bis 15, 

bei der der Prozessor derart eingerichtet ist, daB bei der 
Optimierung der Parameter der Entmischmatrix wenigstens ein 
Diagonalparameter der Entmischmatrix auf einen vorgegebenen 
Wert gesetzt wird. 

17. Anordnung nach einem der AnsprUche 13 bis 16, 

bei der der Prozessor derart eingerichtet ist, daB die 
Entmischmatrix auf eine finite Impulsantwort begrenzt wird. 

18. Anordnung nach einem der AnsprUche 13 bis 17, 

bei der der Prozessor derart eingerichtet ist, dafi die 
Entmischmatrix wShrend der Kumulantenminimierung (5) durch 
Projektion in einen Einheitskreis stabilisiert wird. 

19. Anordnung nach einem der AnsprUche 11 bis 18, eingesetzt 
zur Trennung uberlagerter, statistisch voneinander unabh^ngiger 
Eingangssignale. 
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20. Anordnung nach einem der AnsprQche 11 bis 18, eingesetzt 
zur Trennung tiberlagerter, statistisch voneinander unabhangiger 
akustischer Eingangssignale. 
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Zus aannenf as sung 

Verfahren und Anordnung zur Ermittlung von Par ame tern eines 
technxschen Systems 

Es werden Parameter eines technischen Systems , mit dem 
Ausgangssignale aus einer Menge xiberlagerter, statistisch 
voneinander unabhSngiger Eingangssignale ermittelt werden 
kOnnen, bestimmt. Die Parameter werden derart ermittelt, daB 
die statistische UnabhSngigkeit der Ausgangssignale maximiert 
wird. 
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